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Abstract We study transport through a strongly correlated quantum dot and show that Coulomb blockade 
can appear even in the presence of perfect contacts. This conclusion arises from numerical calculations of 
the conductance for a microscopic model of spinless fermions in an interacting chain connected to each lead 
via a completely open channel. The dependence of the conductance on the gate voltage shows well defined 
Coulomb blockade peaks which are sharpened as the interaction strength is increased. Our numerics is 
based on the embedding method and the DMRG algorithm. We explain the emergence of Coulomb blockade 
with perfect contacts by a reduction of the effective coupling matrix elements between many-body states 
corresponding to successive particle numbers in the interacting region. A perturbative approach, valid in 
the strong interaction limit, yields an analytic expression for the interaction- induced suppression of the 
conductance in the Coulomb blockade regime. 

PACS. 73.23.Hk Coulomb blockade; single-electron tunneling - 73.23.-b Electronic transport in mesoscopic 
systems - 71.27.-|-a Strongly correlated electron systems; heavy fermions 



1 Introduction 

The discreteness of the electron charge blocks the trans- 
port through a small and relatively isolated conductor 
(usually referred as a quantum dot) at sufficiently low 
temperature and bias voltage. This peculiar phenomenon, 
known as Coulomb blockade, is a paradigm in the physics 
of small condensed matter systems. The Coulomb block- 
ade can be lifted by tuning a gate voltage to a point of 
degeneracy, where the energy cost for adding or removing 
an electron to the system vanishes. The obvious signature 
of large conductance oscillations as a function of gate volt- 
age makes this phenomenon promising for applications, 
and at the same time a privileged set-up to investigate 
interactions in confined systems |ll2j . 

A question that has been present since the beginning 
of the studies on Coulomb blockade is: what are the mini- 
mum ingredients for it to be observed? The repulsive inter- 
particle interaction is clearly one of them, as well as the 
condition that the charging energy of the dot be greater 
than the thermal energy fceT or the applied bias eV. Hav- 
ing an almost isolated dot, connected to conducting leads 
by weakly transmitting contacts like tunnel-barriers, has 
been thought to be another essential ingredient. Various 
studies have been undertaken to follow the evolution of the 
Coulomb blockade oscillations as the dot becomes better 
coupled to the leads. 

It was in the quantum Hall regime that Coulomb 
blockade effects in open (well coupled) dots were initially 
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investigated |3I4) . The existence of spatially separated 
edge channels made it possible to observe conductance 
oscillations even in the case in which the point contacts 
at the entrance of the dot had a conductance Gb larger 
than the quantum conductance Go = e'^/h. 

In the absence of the quantum Hall effect, the regime 
of weak Coulomb blockade, which is intermediate between 
that of strong Coulomb blockade (Gb ^ Go) and the one 
of a large number of fully open channels (Gb 3> Go, where 
no Coulomb blockade can be observed), has been studied 
|5lf)| . In the regime of weak Coulomb blockade (Gb ~ Go), 
the precise behaviour of the system depends on the way in 
which the value of the entrance conductance is obtained. 
In the case of metallic islands coupled to leads via tunnel- 
barriers, Gb can be of the order of Go due to a large 
number of weakly transmitting channels; while in semi- 
conducting dots coupled to leads via point contacts, con- 
ductances Gb of the order of Go are obtained by increasing 
the transmission of only a few channels. 

Nazarov has shown that charge quantization may per- 
sist when the dot is connected to the leads via arbi- 
trary conductors |7j, as for example disordered metallic 
wires, instead of tunnel-barriers. In the general case of 
a number of not perfectly transmitting channels, the ef- 
fective charging energy is exponentially small (scaling as 
exp (— Gb/Gq)). However, the charging energy vanishes if 
one of the channels is perfectly transmitting. 

It is important to make the distinction between perfect 
transmission and perfect contact. The notion of transmis- 
sion of the barrier, which should be "measured" between 
two electrodes, may not be a relevant concept when the 



2 



Gabriel Vasseur et al.: Coulomb blockade without potential barriers 



barrier is between a dot and an electrode. In this case 
a perfect contact does not guarantee the absence of re- 
flection, since electrons can be coherently backscattered 
inside the dot. 

Recent works |8I6| considering open dots have taken 
into account the coherent backscattering induced by the 
dot itself. They found that interactions do not modify the 
average conductance when all the channels are either open 
or closed (case of ideal contacts). This conclusion is in 
agreement with an experiment performed on rather large 
chaotic dots [H]. 

Experimentally. Coulomb blockade physics has been 
observed recently in very small silicon based MOSFETs 
where the tunnel-barriers are not built-in [lOlllj . It is 
likely that the Coulomb blockade features of these samples 
are due either to diffusive leads or to electrostatic poten- 
tial barriers created between the dot and the leads. It is, 
however, of interest to study whether it is really necessary 
to have poor contacts for obtaining the Coulomb blockade. 
For dots of reduced dimensionality and containing a low 
number of electrons, we expect electronic correlations to 
be of crucial importance. Therefore we may ask whether 
Coulomb blockade can arise even with perfect contacts, 
due only to strong correlations. 

In order to address this question, it is important to 
go beyond the usual treatment of the interactions in the 
dot based on the capacitive charging energy. Microscopic 
models of correlated chains attached to semi-infinite leads 
have recently been investigated |12ll3ll4j . For good but 
not perfect contacts, it has been found that Coulomb 
blockade-like features are reinforced by the interactions. 

In this paper, we present a study of a microscopic 
model for a strongly correlated quantum dot and show 
that Coulomb blockade arises even in the presence of per- 
fect contacts. 

The model and the method are introduced in section 
121 The numerical results presented in section |2| show that 
Coulomb blockade appears as the interaction strength in- 
creases, even in the absence of potential barriers between 
the dot and the leads. The features of the peaks and valleys 
are analyzed in detail in section ^ using perturbative ap- 
proaches at weak transmission. In the concluding section, 
we discuss the implications of our findings and possible 
extensions of our model in order to approach the current 
experimental setups. 

2 Embedding method for the conductance of 
a one-dimensional model 

Most of the analytical approaches to Coulomb blockade 
rely on the so-called constant charging model, where cor- 
relations are completely disregarded. The only electron- 
electron interaction term considered is the capacitive 
charging energy of the quantum dot. While various re- 
finements, such as the Random Phase Approximation jl5j . 
have been developed to describe the case of small dots with 
poor screening, treating electronic correlations to all or- 
ders is obviously out of reach in any analytical approach. 
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Figure 1. The central region of an infinitely long, one- 
dimensional chain. The particles experience nearest neighbour 
interaction of strength U on sites 1 to Ls- Below the chain 
are sketched the compensating potential of the interaction (see 
text) and the potential due to the gate voltage Vg. 

Numerically, various schemes for obtaining the addition 
spectrum of a quantum dot have been used |15| . However, 
these calculations of ground state energies cannot capture 
the richness and complexity of a transport problem, espe- 
cially for open quantum dots. 

Even numerically, the transport through strongly cor- 
related systems has remained an intractable problem un- 
til recently. The development of the so-called embedding 
method fUJlSj has been an important advance in this 
difficult problem. Relating the transport properties of a 
quantum dot with the thermodynamic properties of the 
combined system (sample -I- leads, as we explain in the 
sequel), allows to extract the conductance through the 
dot from the numerically accessible ground state energy 
of the combined system. 

The embedding method can readily be used for spin- 
less fcrmions (spin-polarized electrons) in one-dimensional 
chains. While the method has been successfully general- 
ized to spin-half electrons [T2|, the problem of spinless 
fcrmions remains more easily tractable. For these reasons 
we restrict ourselves in this work to one-dimensional spin- 
less fermions. It is expected that the qualitative aspects 
of our results apply also in higher dimensions and once we 
take the electron spin into account. 

We consider a model of spinless fcrmions in a chain, 
as sketched in Fig. The dot corresponds to a region of 
length is in which the particles are interacting. In the rest 
of the chain, representing the leads of the standard experi- 
mental set-up, the particles do not interact. Therefore, far 
away from the dot, a Fermi energy Ep is well-defined. We 
fix the latter to the centre of the band, corresponding to 
half filling. The Hamiltonian of the whole system reads 

H = Hk + Hu + Hg, (1) 

where 

oo 

Hk^- + h.c.) (2) 

— oo 

stands for the kinetic energy of an ideal one-dimensional 
chain. Here, Ci annihilates a particle on site i. Setting the 
hopping amplitudes in iJx to unity defines our energy 
scale. 

The dot is represented by the sites i = 1 to Ls (the 
lattice constant is set to unity). The inter-particle inter- 
action 

Ls-l 

Hu = Uj2in^- l/2)(n,+i - 1/2) (3) 

i=l 
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acts only in this region and is assumed to be restricted 
to nearest neighbours (with strength U and fii = c\c^). 
Hu includes a one-body compensating potential (due to 
the presence of the 1/2's, see Fig. Q which ensures, in 
the absence of gate voltage, particle-hole symmetry even 
at /7 ^0. It can be seen as a positive background that 
prevents the particles from leaking out of the dot region, 
in spite of the interaction. Additional particles can then 
be attracted in the dot by the one-body potential 

is 

HG = -V,J2h, (4) 

i=l 

describing the effect of an applied gate voltage Vg. 

The relative isolation of a dot is usually achieved by 
connecting it to the leads via weak links or tunnel-barriers 
(i.e. by choosing a smaller hopping amplitude or a high on- 
site potential) , but we do not take any of these approaches 
in this work. In contrast, in our case, the contacts between 
the dot and the leads are perfect. The dot is therefore only 
defined by the region where the electron-electron interac- 
tions and the gate voltage are present. 

We use the embedding method which allows to de- 
termine the modulus |t| of the effective transmission am- 
plitude t = |t|e'" of the system at the Fermi energy 
of the electrodes, taking fully into account the effects 
of electronic correlations. This thereby yields the zero- 
temperature, linear-response conductance g = |tp (in 
units of e'^/h). The modulus |i| is obtained from the per- 
sistent current flowing in a ring formed by the interacting 
dot together with an infinitely long, non-interacting lead. 

The relation between the persistent current and the 
conductance is based on two facts, which are both realized 
in the limit of an infinitely long lead. First, the persistent 
current in a ring including a one-body (non- interacting) 
scatterer only depends, in this limit, on the modulus |t| 
of the transmission amplitude of the scatterer. Second, it 
is only in the limit of an infinitely long lead that a ring 
containing an interacting part behaves as a Fermi liquid, 
and that a one-body transmission amplitude t can be at- 
tributed to the interacting region attached to leads \2Q\ . 

The study of closed, finite-size systems thus allows to 
obtain the Landauer conductance of the dot in the setup 
described by the Hamiltonian IpJ, provided that an ex- 
trapolation to infinite size is performed. We are able to 
perform reliable extrapolations for an interacting dot of 
Lg = 6, even near the resonances, where the extrapo- 
lation procedure is most difficult |^|. Therefore, besides 
the results for the physics of Coulomb blockade, this work 
also constitutes a test of the embedding method in a com- 
putationally demanding situation. 

To be specific, we use the phase sensitivity T) = 
(L/2)|i?p — Ea\, where i?p,A denotes the ground state 
energies of the ring for periodic (antiperiodic) boundary 
conditions, and L is the total size of the ring. The extrap- 
olation to L oo is carried out from data calculated for 
different ring sizes. The limiting value Poo leads to the ef- 
fective transmission amplitude of the system at the Fermi 
energy |17I21| 

|<| -sin(Poo). (5) 



Similar approaches using the persistent current have also 
been recently proposed |22l28l24j . 

We take advantage of the efficiency of the Density Ma- 
trix Renormalization Group algorithm (DMRG) |25I26| 
which allows to calculate T) with very high precision, tak- 
ing fully into account the electronic correlations. The ex- 
trapolation can then usually be performed from the many- 
body ground state energies of rings of up to 100 sites. Near 
the resonances, however, the extrapolation can be quite 
difficult and we need to treat rings containing up to 300 
sites, keeping up to 1000 states in the DMRG iterations. 



3 Interaction-induced conductance 
oscillations as a function of gate voltage 

The dimensionless conductance g as a function of the gate 
voltage Vg is shown in Fig. |2 together with the mean 
charge inside the dot region. The three graphs shown 
are for a dot of length Lg = 6, and different interaction 
strengths. 

The non-interacting case {U = 0, Fig. |3i) corre- 
sponds to free particles in a chain with a potential well 
of length Ls and depth Vg. The structures in the con- 
ductance for this case can be understood from one-body 
quantum transport through the well potential, and thus 
t{E-p) = |t|e*" can be readily calculated. The conductance 
is then obtained as g = |tp which, when \Vg\ < 2 and for 
-Bp = 0, can be written as 

itp=(cos-(/(Fg)Ls)+ -^"^(;(^y )'\ (6) 

V sm (/(Vg)) J 

with f{Vg) — arccos (— Vg/2). The result, plotted in Fig. 

is in perfect agreement with the numerical data ob- 
tained from the embedding method. This is a stringent 
test of our method since the non-interacting case does 
not represent a special situation for the DMRG algorithm. 
Furthermore, the transmission phase a can be related to 
the number N of particles inside the dot by the Friedel 
sum rule |27I28| 

iVd = a/^, (7) 

where A^d is the number of particles displaced by the dot. 

gives the number N of particles inside the dot only for 
a system with homogeneous density, such that deviations 
between the analytic result (dotted line in Fig.|3t) and the 
DMRG data appear with increasing gate voltage, when 
part of the displaced particles are in the leads. We also 
calculated TV from a direct diagonalization of the one-body 
Hamiltonian, and find perfect agreement with the DMRG 
data (sec Fig. dashed line). 

As one increases the interaction, a well-defined peak 
structure appears. The electron-electron interaction sup- 
presses the conductance between the peaks which become 
sharper with increasing interaction strength. This ten- 
dency can be observed already at moderate interaction 
strength U = 2 (Fig. I^b), and becomes very pronounced 
at strong interaction U = 10, as can be seen in Fig.[2t. In 
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Figure 2. Dimensionless conductance g (triangles) and mean 
number of particles (circles) inside the dot (of Ls = 6 sites) , 
as a function of the gate voltage Vg without interaction (a) and 
Vg/U for [/ = 2 (b) and [/ = 10 (c), obtained by DMRG calcu- 
lations. In the non-interacting case (a), g is known analytically 
(Eq. ®, solid line). The Friedel sum rule (dotted line) gives 
A'^ approximately (see text) while one-body numerics (dashed 
line) is in perfect agreement with the DMRG points. In (b) 
and (c) lines are guides to the eye. 



a one-body problem with symmetric contacts the conduc- 
tance at resonance reaches unity. In our strongly corre- 
lated problem the limited amount of data points and the 
uncertainty of the extrapolation do not allow to settle the 
precise values at resonance in Fig. 

The emergence of a well-defined peak structure in the 
conductance g is accompanied by the appearance of a 
staircase dependence on Vg of the mean particle number 
N inside the dot region. The plateaus appear at integer 
values of the mean particle number, and the steps in N 



occur precisely at the gate voltage values of the conduc- 
tance peaks. 

Although the dot region is perfectly connected to the 
leads, these structures are clearly due to the Coulomb 
blockade effect, and their features can be identified with 
the properties of the usual Coulomb blockade conductance 
oscillations. It therefore appears that the dot becomes ef- 
fectively decoupled from the leads by the electron-electron 
interaction. The underlying mechanism will be discussed 
in the following sections. In particular, wc will show in 
section 14.21 that when the interaction is strong enough, 
the peak positions can be predicted precisely from the 
knowledge of the ground state energies of the isolated dot 
system for different mean particle numbers. Similar con- 
clusions have been drawn from numerical studies of a chain 
containing 7 sites, but connected to the leads by reduced 
hopping matrix elements |13I14| . In this DMRG 
evaluation of the Kubo formula and a DMRG study of the 
real-time dynamics have been performed showing a rein- 
forcement of Coulomb blockade features when the inter- 
action strength is increased. Including the spin degree of 
freedom leads to the appearance of conductance plateaus 
in the Kondo regime. A numerical renormalization group 
(NRG) study [T2| of a Hubbard chain of 4 sites with mod- 
erate interaction strength and reduced coupling to the 
leads shows that the corresponding features arc smoothed 
when the coupling to the leads is increased. 

Since the full Hamiltonian is symmetric under the ex- 
change of particles and holes with a simultaneous reversal 
of the sign of the gate voltage, the results for negative gate 
voltage can be obtained directly from the data presented 
in Fig. 121 The conductance is symmetric with respect to 
the gate voltage, gi~Vg) ~ 9{Vg), while the number of 
particles exhibits the property N{—Vg) ~ Ls — N{Vg). 
Thus, as expected for a dot with Ls = 6 sites, there is 
a total of 6 peaks with the particle number increasing in 
steps from to 6. Systems with larger Lg are therefore 
expected to give qualitatively similar results with Lo, con- 
ductance peaks and Ls charge steps. In the case of odd 
Ls, the number of conductance peaks is odd. Because the 
conductance is an even function of the gate voltage, one of 
the conductance peaks must be at = 0, in contrast to 
the case of even Ls, where = is always in the centre 
of a valley between two conductance peaks. This observa- 
tion allows to understand the even-odd oscillations in the 
conductance of an interacting chain without gate voltage 
[T7| discussed in detail in Ref. [^Hl ■ The oscillations of the 
conductance with the length of the interacting region are 
a very general feature of interacting chains. A perturba- 
tive treatment of the interactions [201 indicates that they 
appear in Hubbard chains as well. They may be related 
to experimentally observed conductance oscillations with 
the length of mono-atomic chains |31| . 



4 Perturbative approaches for weak 
transmission 

In order to develop our understanding of the physical fea- 
tures found numerically, we present an analytical study of 
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the model for the regime of strong interaction. The be- 
haviour observed in Fig. [21 and discussed in the previous 
section, indicates that the dot becomes decoupled from 
the leads when the interaction strength is increased. We 
will show that this is indeed the case. To this end, we treat 
the conductance through the system with a Fermi golden 
rule approach, valid when the effective coupling between 
the dot and the leads is small. In such a framework, the 
starting point is a dot region which is isolated from the 
leads. The coupling between the dot and the leads is then 
treated as a perturbation. In this case, the relevant cou- 
pling is not given by the one-body hopping amplitude be- 
tween the sites 1 and Lg ^ + li which in our 
model is always unity, but by the matrix elements of the 
terms of the Hamiltonian which couple the interacting 
dot region to the semi-infinite leads, sandwiched between 
the full many-body states of the chain. This coupling part 
of the Hamiltonian reads 



He = -(4ci + h.c.) - (cLc^g+i + h.c.) 



(8) 



and corresponds to Hk restricted to the two links between 
the dot and the leads. We now evaluate the lowest order 
contribution in He to the conductance of our system. In 
the case where the involved matrix elements of Hq are 
small (compared to U), the conductance is dominated by 
this term. 

The conductance is related to the transmission of elec- 
trons through the system. Let us start by considering 
Fermi's golden rule for the transition rate 



^\Mf,\'SiEf-E,) 



l^^f ^ —\^^'^fA'^{^f - (9) 

of a particle through the dot region. The initial state 

\i) = Ifci-ft) ® |0(iV)) (10) 

describes a product state built from eigenstates of the iso- 
lated dot and the semi-infinite leads. Here, we take the 
dot in its A^-particle ground state |0(A^)) and an addi- 
tional particle jfc'''^*) with wavenumber k in the left lead. 
In the final state 



I/) = |0(iV)) ® \k 



right \ 



(11) 



the additional particle appears in the right lead while the 
dot returns to its A^-particle ground state. The effective 
transition matrix element is of second order in Hq, and 
can be written as 



'c 



a)(a|gc|/) 



Ei — E, 



(12) 



Here, the sum runs over all intermediate states |a) corre- 
sponding to N+1 particles in the dot region and the Fermi 
vacuum in the leads or to A'^ — 1 particles in the dot and 
both leads occupied by an extra particle. Summing the 
transition rate ^ over the states of the non-interacting 
leads in a small energy interval corresponding to an in- 
finitesimal bias voltage V, one can obtain the correspond- 
ing current and thus the dimensionless conductance 



where p is the density of states in the leads. 

The approach is justified when the effective coupling 
given by the matrix element in Eq. I|12() is small. This is 
the case at large U when the energy denominator is large, 
even though we have perfect contacts and He is not small. 
In addition, the matrix elements of Hq can be reduced due 
to the interaction-induced modification of the dot wave- 
functions. 



4.1 Conductance in the limit of strong interactions 

In the strong interaction limit (U 3> 1,V^), the ground 
state of the isolated dot region for an even dot length ig 
at half filling {N = is/2) is given by a charge density wave 
(or Mott insulator, where the particles occupy alternating 
sites). In this case we can evaluate Mf^i, and hence obtain 
the conductance H13|) in the corresponding conductance 
valley. 

In this section, we treat not only the coupling between 
the dot and the leads, but also all of the hopping terms 
Hk as a perturbation to the other terms Hq = Hu + Hq 
of the Hamiltonian. This corresponds to an expansion in 
1/U and follows the spirit of the large U expansions of the 
persistent current presented in Ref . [IV2\ . In the absence of 
Hk, the two realizations of the Mott insulator with the 
particles on the even sites = CjcJ . . . c^_, |0) and on 

the odd sites \iPq) = cjcg . . . c^g_;^|0) are degenerate. They 
become coupled in A^th order in Hk- The degeneracy is 
therefore lifted even for infinitesimal 1/U and the ground 
state is given by the symmetric (the effective coupling is 
— (1/C/)^^/^) superposition 



V2 



(lV^o) + lV^°))- 



(14) 



(13) 



In the leads, Hk yields plane wave eigenstates with the 
boundary condition that, due to the Mott-Hubbard gap 
(of the order of U), their wave- function vanishes (like 
1/J7) on the first interacting site. For leads of length L 
on either side of the system this gives one-body states 
with energy = —2 cos [k) and wave-functions (z|fc'''^') = 
■\/2/Lsin (fc(l — i)) on the left hand side (i < 1) and 
(i|fc"sht) = y27Isin(fc(i - is)) on the right hand side 
{i > is) of the interacting region. 

Zero-temperature transport through the interacting 
segment of the chain requires that particles are transmit- 
ted elastically. Therefore, we need the lowest order pro- 
cesses in 1/U linking the initial state \i) = |fc'''^*) (8) |^o) 
to the final state |/) = \^po) ® \k"s^*'). These processes 
are the ones in which a particle is transmitted while the 
component |-0q) of the ground state (|14|l of the interacting 
region is connected to the component \iPq) by sequences 
of -I- 1 successive single particle hoppings. In such a pro- 
cess, each of the N particles inside the interacting region 
hops to the next site, and hence the one initially sitting 
on site is leaves the region towards the right lead. In ad- 
dition, a particle from the left lead enters the dot region 
and appears on site 1. An example of such a sequence is 
sketched in Fig. |3| 
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Figure 3. Example of a hopping sequence of Ls/2 + 1 = 4 
hops connecting the two components of the ground state Ijl4|l 
of the uncoupled dot. Full and empty circles represent occupied 
and empty sites, respectively. In this example, a particle enters 
from the left in the second hop, and another one leaves it to 
the right in the last hop. 



These dominating processes in the transport through 
an interacting chain arc similar to the co-tunneling pro- 
cesses through arrays of quantum dots |S] . 

The lowest contributions to the effective matrix ele- 
ment Mf i are therefore of order + 1 in Hk ■ It is given 

by 



S—{ai ,a2 ,. 



{i\HK\ai){ai\HK\a2) 



■ Qjv} 



{Ei — Ea^){Ei — Ea^) ■ ■ ■ (Ei — Ea„) ' 



(15) 

where S = {ai, a2, ■ ■ ■ a^v} are sequences of N intermedi- 
ate eigenstates of Hq linked by iV 4- 1 subsequent single 
particle hoppings. For all the {N + 1)1 different permuta- 
tions of the order of the hops, the numerator is given by 
(— 1)-'^+-'^ sin^ {k) /L. The denominator involves U and Vg 
and depends on the sequence. For small even is we have 
explicitly derived the formula 



1 



sin^ (fc) 



U 



- [/2/4 



is/2 



(16) 



and we have confirmed its validity for larger values of Ls 
by performing the summation over the sequences numeri- 
cally. With H13() and the one-dimensional density of states 
p = L/ {2tt sin (fc)) we get 
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U 



is 



(17) 



for the dominating contribution to the conductance in the 
limit of very strong interaction strength [/, at E-p = 
(fcp = Ti'/Z). This demonstrates that Coulomb blockade 
eventually occurs for sufficiently strong interactions. The 
suppression of the conductance in the valley due to the 
Mott insulating behaviour follows a power law in the in- 
teraction strength U and is exponential in the dot length 
Ls ■ The result of Eq. H17|) is a generalization of the result 
found for = in ^H]. The small values of the conduc- 
tance around Vg = in Fig. ^ arc well described by Eq. 
H17(l and the agreement improves for stronger interaction. 

The power-law decay in Eq. H17|) of the valley con- 
ductance with very large interaction strength U is ob- 
tained analytically and confirmed by very accurate nu- 
merics in our case of spinless fermions. For the case of 
Hubbard chains, an exponential decay was concluded |33| 
from NRG data for relatively moderate values of the inter- 
action. However, our analytical approach can be applied 
to the Hubbard case as well, leading to a power law de- 
pendence of the conductance on the interaction strength 



in the limit of large interaction strength. This power law is 
a consequence of the lattice model. The exponential decay 
mentioned in [33j might be an apparent dependence in a 
regime of intermediate interaction strength. 

Within the same model, but in the absence of a gate 
voltage, it was previously found ^7] that the smoothing 
of the contacts (i.e. the slow branching of the interaction 
going from the leads to the dot) reduces the interaction- 
induced suppression of the conductance found for even 
values of Ls- In the framework of the present paper, this 
means that smoothing the contacts will reduce the depth 
of the conductance valley at zero gate voltage. However, an 
adiabatic branching of the interaction is necessary to re- 
cover the perfect conductance found without interaction. 
Furthermore, numerical results |34j show that smoothing 
the interaction does not necessarily suppress the other val- 
leys. 



4.2 Peak positions and shapes 

From the data presented on Fig. |21 the positions V^^^'' 
of the resonances can be determined as the gate voltage 
values for which there are precisely TV — 1/2 particles inside 
the dot. These peak positions V^^^-*, in units of U, are 
shown in Fig. 0] (symbols) as a function of the interaction 
strength U, together with the positions 



yW(t^) - Eo{N, U)-Eo{N~l,U), 



(18) 



that are obtained from the many-body energies Eq{N, U) 
of the isolated dot with N particles (lines). One can see 
that already for [/ = 10, the positions of the resonances 
are very close to the isolated dot prediction of Eq. H18(l . 
This is an additional confirmation of the fact that the 
interaction leads to a decoupling of the dot from the leads. 

One of the most relevant features of the numerical re- 
sults shown in Fig.[21is the sharpening of the conductance 
peaks with increasing interaction strength. We shall now 
provide an understanding for the effect of strong interac- 
tions on the peak width. 

The width Jjv of a peak is given by the value of the 
matrix element of Hq. in the numerator of H12(l . Since the 
conductance resonance at Vj,^^-* is given by the degeneracy 
of the N- and N — 1-particle ground states, we restrict our 
discussion to the contribution of the transition between 
those two ground states. Within a qualitative discussion, 
we can characterize the parameter dependence of the peak 
width using the transition matrix element 

-(fc{?f*|®(0(iV-l)|cSci|0(7V)), (19) 
which contains the amplitude 

= (0(iV-l)|ci|0(7V)). (20) 

The main contribution to Pjv is proportional to \'^ion\^ ■ 
This matrix element depends on the many-body wave- 
functions, and is therefore strongly influenced by the in- 
teractions. For general values of A'', it can be highly non- 
trivial to determine the ground state of the isolated dot. 
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Figure 4. Positions Vg^^' /U of the peaks from the numeri- 
cally obtained conductance (as in Fig. |5J as a function of the 
interaction strength U (symbols). The lines give the isolated 
dot prediction V^^'^ /U according to Eq. H18|l . The inset shows 
the peak width Fm as a function of U in double logarithmic 
scale. The error bars are given by the uncertainty of the fitting. 
For TV = 4 they are smaller than the symbol size. The A'^ = 6 
data points (not shown) are the lowest and exhibit the largest 
error bars. The dashed line represents the asymptotic power 
law 

There are, however, particular values of TV for which we 
can study analytically. 

One of these special values in the case of even Ls is 
TV = Ls/2 + 1, which corresponds to the left peak in Fig. 
13 In the limit of strong interactions, the |0(TV — 1)) state 
is a superposition of two charge density wave configura- 
tions, as presented in Eq. H14|l . The additional particle 
in |0(TV)) is a defect (with interaction energy U) of the 
Mott insulator which can propagate freely, and can be 
described by an eflfective one-particle theory. In the ab- 
sence of the kinetic energy ffx, the state |0(TV)) is a su- 
perposition of the Ls /2 degenerate low- lying states of the 
form cjcg . . . c\^_■^c\^ . . . c^^ |0). In the presence of Hk, the 
weights of the different components are obtained by solv- 
ing the equivalent problem of a free particle on Ls/2 sites. 
Only the weight of the first state (namely for i = 1) enters 
the expression for '^Ls/2+Ij one obtains 

In this limit of strong interaction, the lowest order contri- 
bution to '^Ls/2+i is thus independent of U (the limiting 
value is reached already for [/ > 2 as can be seen in the 
inset of Fig.0J, but decreases as Lg ^^'^ with the length of 
the system when Ls is large, indicating a decrease of the 
peak width proportional to Lg ^. 

It is of some interest to compare this behaviour at 
large Lg with a non-interacting evaluation of the matrix 
element which one can calculate although at ?7 = the 
dot is not decoupled from the leads. In this non-interacting 
case, the eigenstates are given by free particles on Lg sites. 
^Ls/2+i is then obtained from the value of the (Ls/2+l)th 
lowest one-body wave-function on the first site. This leads 

— 1/2 

to a size dependence of the matrix clement as Lg only. 



From this, one may conclude that the interactions play an 
important role in narrowing the resonance. 

The other case in which analytical calculations are 
straightforward is the one of the last peak when TV = Lg . 
In this case |0(TV)) corresponds to the full dot. The state 
|0(TV— 1)) corresponds to one hole in a potential well with 
U/2 steps on the edges (as sketched in Fig. It is ob- 
vious that when U increases, the hole is less likely to be 
found on an edge site, because of the potential step. The 
ground state of such a one-body system is solved by the 
ansatz of an even wave- function with '0(l<n<Ls) = 
Acos(fc(n- (Ls + l)/2)) and = ^{Ls) = B. The 

normalization allows to get B, which yields 




where k is the smallest positive value that satisfies the 
condition 

cos " ~T' ^'^^ (~V~"^) ■ ^'^^^ 

In this case, it can be seen that ^ depends on U via k. In 
the limit of strong interaction, 1)23(1 can only be fulfilled 
if cos (■^^^— '-fc) ~ U^^ . This implies that the numerator of 
((22|l decreases proportionally to and thus the peak 
width Fn is expected to decrease as U~^. This impressive 
example of the reduction of the peak width due to the 
interaction is illustrated in the inset of the Fig. 0] 

We end this section with a discussion on the symme- 
try of the peaks. As can be seen in Fig. Et, the conduc- 
tance peaks are not of symmetric Lorentziaii shape. This 
is most obvious for the first peak and can be explained by 
an asymmetry in the coupling of the many-body states to 
higher and lower mean particle numbers. In order to un- 
derstand this, we consider the representation of the matrix 
element Mf^i as a single sum over the eigenstates \a) of 
the isolated dot, as in Eq. ((T^ . 

Outside resonances, the number TV of particles inside 
the dot, and therefore the initial (|10|) and final states, 
are well-defined, and one can use Eq. H13() to evaluate the 
conductance. Very close to the conductance peak however, 
the transmission becomes large and the Fermi golden rule 
approach is not valid. To address the issue of the asym- 
metry of the peak, we consider values of the gate voltage 
which are in the vicinity of a resonance V^^^^ (and away 
from all the other) , but far enough from the peak centre to 
have a small conductance. In this way, the conductance in 
the tails of the peaks can be treated using ((12|l and H13() . 

The main contributions to the sum (|12|l depend on 
whether the gate voltage is above or below the resonance 
value Vg^^\ This resonance corresponds to the transition 
from TV — 1 to TV particles in the dot. 

For gate voltage values below the resonance, the states 
|i) and |/) are TV — 1-particlc states, and the most impor- 
tant contributions to Mf^i come from the sum over the 
states I a) with TV particles inside the dot. The terms with 
TV — 2 particle states \a) are suppressed by the much larger 
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energy denominator in H12(l . For gate voltages slightly 
above V^''^'', the initial and final states are A^-particle 
states and the main contributions to the sum are due to 
intermediate N ~ 1-particle states. 

Since the many-body excitation spectrum of the iso- 
lated dot can depend dramatically on N, the sum can give 
very different results on the two sides of the peaks, which 
are therefore expected not to be symmetric. 

We illustrate this mechanism for the example of the 
first peak of Fig. [2t, corresponding to the transition be- 
tween half fiUing, - 1 = Ls/2 = 3, and = 4. On the 
left-hand side of this peak, the sum is dominated by the in- 
termediate 4-particle eigenstates of the isolated dot, while 
the evaluation of the conductance on the right-hand side of 
the peak is dominated by the 3-particle eigenstates of the 
dot. The 3-particle spectrum of the dot is characterized by 
the Mott-gap which is of the order of U. The contribution 
of the excited 3-particle states is therefore suppressed by 
energy denominators of the order of U, leading to a smaller 
value of Mf^i on the right-hand side of the peak than on 
the left-hand side, thus explaining the peak asymmetry 
observed on Fig. Et. 



5 Conclusion 

We have shown that Coulomb blockade physics can occur 
when an interacting system is coupled to leads by perfect 
contacts. The underlying mechanism is that while the con- 
tacts are perfect on the one-particle level, the interaction 
can introduce many-body effects which effectively reduce 
the coupling between the interacting system and the leads. 

In order to illustrate this peculiar behaviour, we 
have studied the gate-voltage dependence of the conduc- 
tance through a one-dimensional chain in which spinless 
fermions are strongly interacting. We have used the em- 
bedding method to compute the conductance g and the 
mean number N of electrons inside the interacting chain 
from numerical DMRG data for the many-body ground 
state properties of large rings. This approach takes into 
account all the electronic correlations. Despite the absence 
of potential barriers separating the interacting chain from 
the leads, the numerical results in the regime of strong 
interaction unambiguously show features which are char- 
acteristic of Coulomb blockade physics in the transport 
through quantum dots. Peaks appear in the gate-voltage 
dependence of which are accompanied by steps in N. 
These structures become more and more pronounced, and 
the conductance becomes more and more suppressed in 
the valleys between the peaks, as the interaction strength 
is increased. 

Our results indicate that, in spite of the perfect con- 
tacts, many-body effects lead to an effective decoupling 
of the interacting chain from the leads when the inter- 
action is sufficiently strong. One can then consider the 
interacting part of the chain as a quantum dot. In or- 
der to interpret these results, we presented an analytical 
study of the model. We explained the deepening of the 
valleys (and therefore the occurrence of Coulomb blockade 



itself) from a perturbative calculation of the conductance 
in the regime of strong interaction. In addition, we showed 
that the positions of the peaks can be deduced from the 
many-body eigenenergies of the isolated dot at strong in- 
teraction, confirming that the dot becomes effectively de- 
coupled from the leads. Furthermore, we explained the 
narrowing of the width of the peaks with the interaction 
strength, as well as the asymmetry of the peak shape, by 
analyzing the effective matrix elements for the dominating 
transitions. 

The current tendency towards the development of 
smaller MOSFETs is likely to lead to regimes where 

the theoretical concepts discussed in this work could be 
properly applied. It would be interesting to further de- 
velop our model in order to bridge the gap with those 
experiments. Addressing the effects of finite temperature, 
longer range interactions, and electron spin are envisioned, 
as well as going beyond the one-dimensional case. The 
newly developed method of obtaining non-linear conduc- 
tances from time-dependent DMRG 14, could be applied 
to our model, and it would be interesting to study the 
suppression of Coulomb blockade by the bias voltage. 

We acknowledge useful discussions with M. Hofheinz, G.-L. 
Ingold, R. A. Molina, J.-L. Pichard and M. Sanquer. We thank 
P. Schmitteckert for discussions and for providing his DMRG 
code. 
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